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Table 1: Details of the unquenched lattice simulations which have been used for the studies of kaon physics. 
The relevant kaon observables computed in these calculations are listed in the last column. 

1. Introduction 

Kaon physics has always played a crucial role for our understanding of fundamental inter- 
actions in the flavour sector. Together with B-physics, for which however precise experimental 
information has become available only in the last decade, the study of kaon physics has allowed 
fundamental tests of the Standard Model and provides, at present, precise constraints on its possible 
new physics extensions. 

While there are few selected processes in kaon physics, like the rare K — > nvv decays, which 
can be studied with almost negligible theoretical uncertainties, in most of the cases the extraction 
of the physical results also relies on our capability of controlling the non-perturbative effects of the 
strong interactions and, therefore, on the accuracy of lattice QCD calculations. This is the case, for 
example, of the determination of the CKM matrix element V IIS from the study of semileptonic and 
leptonic kaon decays, or the theoretical prediction of the Ek parameter which controls the amount 
of indirect CP violation in IC* — K° mixing. The hadronic parameters entering these processes, 
namely the vector form factor /+(0) of semileptonic kaon decays, the leptonic decay constant fx 
and the kaon bag parameter Bk, are the quantities to be computed in lattice QCD calculations. 

In the last few years, the accuracy of the lattice predictions for kaon physics observables is 
significantly increased. Extensive unquenched lattice simulations have been performed by various 
lattice collaborations, using different lattice methods (i.e. different actions, renormalization tech- 
niques, etc.). A list of these simulations, their main details and a compilation of kaon observables 
which have been studied with them, is presented in table [I]. These simulations typically involve 
ensembles with different lattice spacings, which allow the continuum extrapolation to be eventu- 
ally performed. They also include a number of relatively light simulated quark masses, with the 
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lightest pion masses now typically smaller than 300 MeV. In this mass region, a controlled chiral 
extrapolation to the physical light quark masses can be performed, using in most of the cases chiral 
perturbation theory (ChPT) as a useful guideline for the extrapolation. 

The relative abundance of unquenched lattice calculations, in which most (if not all) the sys- 
tematic uncertainties are kept well under control and which use different approaches, characterized 
by different systematics, has allowed to reach in the determination of the corresponding hadronic 
parameters a significantly improved accuracy. It should be also noted that a table similar to table |] 
for lattice studies of B-physics would present, today, a much shorter list. 

In this talk I will mainly concentrate on reviewing the recent lattice results for the three 
hadronic parameters which are of particular interest for kaon physics, namely the vector form 
factor /+(0) of semileptonic kaon decays, the ratio fjclfx of leptonic decay constants and the kaon 
bag parameter Bk- I will also mostly rely, for this task, on the detailed work which is being per- 
formed by the FLAG group, that will be introduced in the next section. For each lattice calculation, 
a colour code rating in the FLAG style will be assigned, and I will also present my best averages 
for the hadronic parameters. I will conclude this talk by briefly summarizing the status of lattice 
studies of non-leptonic K — > Kit decays. 

2. The FLAG working group: rating and averaging criteria for lattice results 

The improved control of systematic uncertainties achieved in the last few years by lattice QCD 
calculations, particularly for observables in the kaon sector, simplifies the task of deriving the 
corresponding lattice averages to be used in phenomenological analysis. This kind of task is one of 
those which is currently being addressed by the Flavianet Lattice Averaging Group (FLAG) [jl]], a 
working group of the Flavianet European network constituted in November 2007. 

The aim of FLAG is to provide for each considered quantity, to the network's working groups 
and to the wider community, the following information: i) a collection of current lattice results and 
references; ii) a summary of the essential aspects of each calculation; iii) averages of lattice results. 

The quantities which are being considered in the first FLAG report [jl]] are the light quark 
masses (m u , m^, m s ), the SU(2) and SU(3) low energy constants, the kaon semileptonic form factor 
/+(0), the ratio of leptonic decay constants fx/ fit and the bag parameter Bk- It is clear from this 
list that a significant overlap exists between the FLAG work and the task I have been given at this 
conference. I will then take advantage of this overlap, and I will use for this talk several FLAG 
results. I am indebted and grateful for that to my colleagues in the FLAG group. 

One of the FLAG proposals that I'm going to follow in this review concerns the way of sum- 
marizing the essential aspects of each lattice calculation. This is done by using an easy-to-read 
"colour code" classification. Specifically, a number of sources of systematic errors are identified 
and a colour with respect to each of these is assigned to each calculation. The prescription is as 
follows: 

*k when the systematic error has been estimated in a satisfactory manner and convincingly 

shown to be under control; 
• when a reasonable attempt at estimating the systematic error has been made, although this 

could be improved; 

■ when no or a clearly unsatisfactory attempt at estimating the systematic error has been made. 
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It should be clear that the precise criteria used in determining the colour coding are unavoid- 
ably time-dependent. The sources of systematic error and definitions which are currently adopted 
by FLAG, and which I'm also going to follow in this talk, are: 

• Chiral extrapolation: 

* M^ min < 250 MeV 

• 250 MeV < M nMn < 400 MeV 

■ M nMn > 400 MeV 

It is assumed that the chiral extrapolation is done with at least a three-point analysis. In case 
of nondegeneracies among the different pion states M n stands for an average pion mass. 

• Continuum extrapolation: 

* 3 or more lattice spacings, at least 2 points below 0. 1 fm 
2 or more lattice spacings, at least 1 point below 0. 1 fm 

■ otherwise 

It is assumed that the action is 0(a)-improved. The colour coding criteria for non-improved 
actions change as follows: one lattice spacing more needed. 

• Finite- volume effects : 

ir {M n L) m [ n > 4 or at least 3 volumes 

• {M n L) inj n > 3 and at least 2 volumes 

■ otherwise 

It is assumed that L m j n > 2 fm, otherwise a red dot will be assigned. In case of nondegenera- 
cies among the different pion states M n stands for an average pion mass. 

• Renormalization (where applicable): 

* non-perturbative 

2-loop perturbation theory (with a converging series) 

■ otherwise 

• Running (where applicable): 

non-perturbative 

• otherwise 

■ — 

Of course any colour coding has to be treated with caution and it goes without saying that 
these criteria are subjective and evolving. Moreover, the extent to which each source of systematic 
uncertainty affects the lattice calculation is observable dependent. FLAG believes, however, that 
this attempt to introduce quality measures for lattice results, in spite of being necessarily schematic, 
will prove to be a useful guide. 

The other main purpose of FLAG is to provide averages of lattice results. The average should 
only include, as far as possible, only "good quality" lattice calculations. This is implemented, in 
practice, by relying on the colour coding: unless special reasons are given for making an exception, 
the averages are restricted to data for which the colour code does not contain any red dot. In 
deriving the averages quoted in this talk, I will follow the same criterium. 

There are two other criteria adopted by FLAG for computing the averages which, however, 
I'm not going to apply for the purposes of the present review. One is related to the publication 
status, for which a coloured symbol is also introduced: 
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• Publication status : 
A published or plain update of published results 

preprint 
C conference contribution 

The FLAG policy is to consider in the averages only calculation which have been published, i.e. 
which have been endorsed by a referee. While I find this policy perfectly justified for the FLAG 
purposes, I also believe that the same criterium is not suitable for the reviewer at the lattice confer- 
ence. The latter is asked to concentrate the attention mainly on the new (and typically unpublished) 
results presented at the conference. It is the task of the reviewer, rather than of an external referee 
in this case, to judge the quality and the reliability of the presented results. 

The other criterium adopted by FLAG is related to the number of flavours of dynamical quarks 
introduced in the simulation. The policy that is being followed by FLAG consists in presenting 
separate averages for the Nf = 2 and Nf = 2 + 1 calculations. This issue has been quite debated 
within the working group, and I personally do not consider the choice currently pursued by FLAG 
an optimal one for a review of lattice results. There are mainly three reasons for that, in my opinion, 
that I would like to mention here, since the issue is also relevant for the present review. 

i) I believe that it would be useful to present, to the wider community of particle physics, the 
"best lattice result" in terms of just a single number, rather than two. When separate averages for 
the Nf = 2 and Nf = 2+1 results are quoted, the natural interpretation of the latter as the best result 
may not always correspond to the actual situation. The error due to the quenching of the strange 
quark is rather small in most of the cases, whereas other sources of systematic uncertainty could be 
better under control in the Nf = 2 determination. 

ii) There are cases, like the one of the semileptonic form factor discussed in the next section, 
in which the error due to the quenching of the strange quark in the Nf = 2 calculation is evaluated 
and included in the systematic uncertainty. In other cases, the comparison between the Nf = 2 
and Nf = 2 + 1 results shows, a posteriori, that the systematic effect due to the quenching of the 
strange quark is not visible within the currently reached accuracy (to the best of my knowledge, 
this is actually the case of all lattice calculations performed so far). In all these cases, I do not 
see any valid reason for ignoring the information coming from the Nf = 2 calculations and for not 
combining together the two sets of results. 

iii) The N f = 2+l theory is not really "full QCD". Indeed, N f = 2 + \ + \ lattice calculations 
are already being performed (see for instance [Q]). I would find unreasonable to simply forget the 
Nf = 2 + 1 calculations when Nf ■ = 2 + 1 + 1 results will be available. For all those quantities for 
which the error due to the quenching of the charm quark will turn out (a posteriori) to be negligible, 
I will suggest again to average together Nf = 2+1 and Nf = 2+ I + I results. 

For these reasons, in deriving lattice averages for the present review, I'm not going to follow 
the FLAG criterium as far as the number of dynamical flavours is concerned, and I will rather 
address this issue on a case by case basis. 

3. |V US | from semileptonic and leptonic kaon decays 

The most accurate determinations of the Cabibbo angle, or equivalently the CKM matrix ele- 
ment Vus, come from the study of semileptonic K — > %tv (K(j) and leptonic K — > £v (K^) decays. 
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Very precise experimental measurements of the K& and of the ratio of Kgx over Ka decay rates 
allow to determine the following combinations of CKM and hadronic parameters 



\Vus I /+(0) =0.21664(48) 



V w 



v, 



ud 



^=0.27599(59), (3.1) 
f% 



with an accuracy of about 2%c. In eq. ([O]), /+(0) is by convention the form factor for the K° — > 
7i~ matrix element. Moreover, both /+(0) and the ratio of decay constants fx/f% in eq. (3.1) 
are defined in the isospin symmetric limit m u = (keeping the kaon and pion masses to their 
physical value) and neglecting electromagnetic corrections. These are therefore the quantities that 
are directly determined in lattice QCD simulations. 

A determination of the hadronic parameters /+(0) and fx/fn, which assumes the validity 
of the Standard Model and it is independent of lattice QCD calculations, has been provided by 
FLAG [p]]. It makes use of the first row unitarity constraint on the CKM matrix, 

\V ud \ 2 + \V us \ 2 + \V uh \ 2 = \ , (3.2) 

which holds in the Standard Model. Since within present uncertainties the contribution of \V u b\ 



in eq. ( p.2| ) is numerically negligible, the unitarity constraint ( p.2| ) can be combined with the two 



experimental results in (3T) to provide a set of three equations and four unknowns: |V„d|, \Vus\, 
/+ (0) and fx/ fit ■ A precise determination of | V u d | is provided by the study of nuclear j8 -decays [Q] , 
based on the analysis of 20 different superallowed transitions: 

1^1=0.97425(22). (3.3) 
Using this result, one can then obtain a determination of the other three parameters, 

\V u .s\ = (1 - \Vud\ 2 f /2 = 0.22547(95) (3.4) 

/+(0) =0.9608(46) , ^ = 1.1925(56) . (3.5) 

Jit 



The estimates (3.5) of the hadronic parameters are benchmarks for the lattice results reviewed in 
this talk. Since both /+(0) and fx/ fit are equal to 1 in the SU(3)-symmetric limit, what it is actually 
measured on the lattice are the SU(3) breaking effects, i.e. the differences /+(0) — 1 and fx/ fit — 1- 



Eq. ( p.5| ) shows that, in order to provide a significant test of the Standard Model, these differences 



must be determined on the lattice with an accuracy of about 10% and 3% respectively. Of course, 



when exploiting new physics scenarios beyond the Standard Model, the unitarity relation (3.2) 



should not be assumed and the determinations in eq ( |3.5[ ) are no longer valid. 

3.1 Semileptonic kaon decays: f+(0) 

The lattice determinations of /+(0) and, in the next section, of fx/ fit, are now reviewed. 

The standard approach to study the vector form factor of decays is based on SU(3) ChPT 
In this framework, the vector form factor at zero momentum transfer has an expansion of the 
form /+(0) = 1 +fi +/4 + • . ., where /„ = 0[m" K J (4nf K )"] and the first term of the expansion 
is equal to 1 due to the current conservation in the SU(3)-symmetric limit. The Ademollo-Gatto 
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theorem ^ shows that the deviation of /+(0) from 1 is at least quadratic in the breaking of SU(3). 
Moreover, the first correction fa receives contribution only from chiral loops and can be computed 
unambiguously in terms of the kaon and pion masses and the pion decay constant. It takes the 
value fa = —0.0226 The problem of estimating /+(0) can be thus re-expressed as the problem 
of finding a prediction for A/ = /+(0) — (1 + fa)- The reference estimate for this quantity is still 
the one obtained by Leutwyler and Roos in 1984 [f7|], using a general parameterization of the SU(3) 
breaking structure of the pseudoscalar meson wave functions. It reads Af = —0.016(8), which 
implies /+(0) =0.961(8). 

Lattice QCD studies of Kb decays started only relatively recently. The strategy, which allows 
to reach the required percent accuracy in the determination of the vector form factor, has been 
developed in ref. [|8p, and it is based on the calculation of the scalar form factor at maximum 
momentum transfer {q 2 nax = (niK — m^) 2 ) through the ratio 

(n\sybu\K)(n\sybu\K) (m K + m n ) 2 2 2 2 
{7i\uYou\7i){K\sYos\K) 4m K m n Ju ^ mwc: n > 



where the external states are of pion and kaon at rest. The double ratio on the l.h.s. of eq. ( p.6| ) 
is equal to 1 in the SU(3) symmetric limit, and it can be evaluated on the lattice with sub-percent 
accuracy at the simulated values of pion and kaon masses. The physical form factor at zero mo- 
mentum transfer, /+(0) = /o(0), is then obtained by extrapolating fa{q^n ax , m 2 K , m 2 ) to q 2 = and 
to the physical meson masses. 

A list of lattice results for / + (0) is collected in table ||, with the colour code in the FLAG style 
assigned for each calculation. The relevant simulation parameters for these calculations are given 
in table ||} The first lattice calculation of the form factor, by the SPQ CD R collaboration was 
performed in the quenched approximation and with rather large values of simulated pion masses, 
m K > 500 MeV. It gave the result /+(0) = 0.960(5)(7), in remarkable agreement with the quark 
model prediction by Leutwyler and Roos [f7|]. Subsequently, unquenched calculations [11 , 12|, jOy , 



performed with Nf = 2 dynamical fermions but still large values of simulated pion masses, con- 
firmed the quenched result. The first lattice calculation of / + (0) which aims to keep under control 
all sources of systematic uncertainties has been performed by the RBC/UKQCD collaboration [S]. 



This year, a second calculation of comparable accuracy has been presented by ETMC [10]. Both 
the update of the RBC/UKQCD 07 result and the new ETMC 09 calculation have been discussed 
in the parallel session at this conference [[14], 15]. 



The RBC/UKQCD calculation [Q] uses the DWF action with Nf = 2+1 dynamical quarks, and 
pion masses as light as 330 MeV. The colour code rating for this calculation, displayed in table ||, 
includes a red square for the continuum extrapolation since the simulation has been performed at a 
single value of the lattice spacing (a ~ 0.11 fm). A 4% of systematic error on (1 — / + ) due to dis- 
cretization effects is assigned on the basis of the parametric estimate Odiscr. ~ (^qcd) 2 , assuming 



Aqcd - 300 MeV. In the parallel talk at this conference |]14|], James Zanotti for RBC/UKQCD has 
announced that a simulation at a second finer lattice spacing (a ~ 0.09 fm) is being performed and 
corresponding results for the vector form factor will be presented soon. Two other improvements 
on the existing calculation are being implemented, namely a simulation at a second value of the 
strange quark mass (the value of m s firstly simulated by RBC/UKQCD [Q] was about 15% heavier 
than the physical strange quark) and the use of twisted boundary conditions. The latters, applied to 
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Ref. 
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/ 
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RBC/UKQCD 07 I ||] 2+1 A 



0.9644(33)(34)(14) 
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QCDSF 07 
RBC 06 
JLQCD 05 



@ 2 A • • 

© 2 C ■ * 

@ 2 A ■ * 

H 2 c ■ * 



0.9560(57)(62) 
0.9647(15)** 
0.968(9)(6) 
0.967(6) 



SPQ CD R04 



0.960(5)(7) 



Table 2: Colour code rating in the FLAG style and summary of lattice results for /+(0). The results dis- 
cussed at this conference are framed. 



Collaboration Ref. Nf action 
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Table 3: Parameters of the simulations listed in table ||. 



valence quark fields [16, [17|, 18], allow to simulate very close or directly at q 2 = 0, thus removing 
the model dependence associated with the momentum extrapolation. 

The ETMC calculation JTo| ] of /+(0) has been presented at this conference by Silvano Simu- 



la Q15|]. It is based on the set of Nf = 2 simulations performed by ETMC with twisted mass fermions 
at maximal twist. Finite size effects on the form factor have been estimated by simulating on two 
different volumes, whereas discretization errors have been evaluated by performing, for a single 
value of the light quark mass, calculations at three different lattice spacing (a ~ 0.07, 0.09, 0.10 
fm). 1 In order to get results close to q 2 = 0, twisted boundary conditions on the valence quark 
fields have been implemented. Both pole and quadratic fits have been considered to interpolate to 



The results on the coarsest lattice have been analysed after the conference and are discussed in the proceedings [151] 
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q 2 = 0. In the main simulation on the lattice with a ~ 0.09 fm, 6 different values of the light quark 
masses have been simulated, with the lightest pion mass being m n ~ 260 MeV. 

The relatively large number of light quark masses simulated by ETMC allows to achieve a 
good control over the chiral extrapolation, which represents one of the main source of systematic 
uncertainty in the lattice evaluation of /+(0). The chiral extrapolation has been performed using 
both SU(3) and, for the first time, SU(2) ChPT. The use of SU(2) ChPT applied to kaon observa- 



bles [19] has been mainly advocated, in the context of lattice calculations, by RBC/UKQCD [EG]. 



For a general discussion on the applicability of the different versions of ChPT to lattice results see 



the plenary talk by Enno Scholz at this conference [ 22 ] . For the vector form factor /+ (0) , the SU(2) 
chiral expansion has been derived at the NLO by Flynn and Sachrajda pi]]. ETMC finds that the 
convergence of the SU(2) expansion for / + (0) is indeed very good, and contributions beyond NLO 
are small both at the physical point and in the region of simulated lattice data. Instead, both ETMC 
and RBC/UKQCD observe that the contribution beyond the NLO in the SU(3) chiral expansion, 
i.e. the quantity Af = /+(0) — (1 + fa), is large and of the same size of the NLO contribution fa. 
Nevertheless, the analyses based on SU(2) and SU(3) ChPT performed by ETMC leads eventually 
to consistent results, and their difference is included in the final estimate of the systematic error. 

The ETMC calculation is performed with Nf = 2 dynamical flavours. The effects due to the 
strange quark loops are exactly accounted for, in the calculation, up to the NLO in the SU(3) chiral 
expansion. This has been possible, because the NLO term in the chiral expansion, fa, can be 
precisely computed (in terms of pion and kaon masses only) for the theories with Nf = [0], Nf = 



2 [23] and Nf = 2+1 [g] dynamical quarks. Thus, the only uncertainty due to the quenching of 
the strange quark in the ETMC calculation concerns the NNLO contribution Af. For this quantity, 
the relative uncertainty has been estimated by ETMC to be of the order of 13% which turns out to 
be of the same size of the entire difference (Af) N f =1 — (Af) N f =0 , evaluated using the quenched 
estimate of Af of ref. [§§]. This difference is most likely an overestimate of the true error affecting 
Af in the Nf = 2 calculation, i.e. (A/)^- - {Af ) N f= 2 . 

In order to quote a lattice average for /+(0), I will take into account both the RBC/UKQCD 
and ETMC results, i.e. 

/+(0) = 0.964(3) (4) [Nf = 2 + 1 , RBC/UKQCD 07 ] 
/+(0) = 0.956(6) (6) [N f = 2, ETMC 09 ] . 

These are the only results which are obtained with relatively light pion masses and do not get 
therefore in table || a red square for the chiral extrapolation. The RBC/UKQCD calculation has a 
red square assigned for the continuum extrapolation, but this source of error is subdominant in the 
calculation (even if the estimated error of 4% were doubled, the impact on the final result would 
be small). On the other hand, except for the partial quenching, the Nf = 2 calculation by ETMC 
has at present a better control over the other systematic uncertainties with respect to the Nf = 2+1 
calculation by RBC/UKQCD (3 lattice spacings rather than 1, 6 pion masses rather than 4, both 
SU(3) and SU(2) chiral extrapolations). Since the uncertainty due to the quenching of the strange 
quark has been accounted for in the ETMC result, I quote as the best estimate of / + (0) the average 



of the two results in eq. ([3/7]), obtaining 

(3.8) 



/ + (0) = 0.962(3)(4) 
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Figure 1: Summary of results for the vector form factor /+(0) ofkaon semileptonic decays. Lattice results, 
as obtained from simulations with Nf = 0,2,3 dynamical quarks, are shown in the upper side of the plot by 
red squares. The results of analytical model calculations /f/|?,/[26|/-j[|2S|/ are shown in the lower side of the plot 
by green triangles. The black circle represents the lattice-independent estimate of /+(0) given in eq. 
which is only valid in the Standard Model. The blue square and the grey band show the average of lattice 



results for /+(0) derived in eq. (31 ). This plot has been produced by FLAG and updated for the lattice 
conference. 



The first error in eq. ( J3 . 8[ ) is statistical, and it is evaluated in the standard way assuming gaussian 
statistical uncertainties. The second error is systematic, and it has been taken to be equal to the 
systematic error quoted by RBC/UKQCD. 2 



The lattice average (3.8) is in very good agreement with the lattice-independent estimate 



quoted in eq. (3.5), based on CKM unitarity and the determination of V U( j from nuclear /3-decays. 



Combining the lattice determination (3.8) with the experimental result of eq. (3.1), one obtains the 
estimate 

\V m \ Kn =0.2252(13) , (3.10) 



in good agreement with the unitarity determination in eq. ( |3.4| ). 

A summary of lattice results for / + (0) is shown in fig. [I], which is an updated version for the 
lattice conference of a plot produced by FLAG [jl|]. The blue square and grey band in the plot 



2 While this contribution was almost finished to be written, an updated estimate of the form factor has been presented 
by RBC/UKQCD Q, which reads 



(3.9) 



/ + (0) = 0.9599(34) (;^)( 14) . 

The new analysis includes the simulation at a second value of the strange quark mass and twisted boundary conditions 
to simulate directly at q 1 = 0, as anticipated at this conference. With respect to ref. [[]], the change in the central value, 
as well as the second error in eq. (3.9 1, are due to the uncertainty in the chiral extrapolation, which was not considered 
in ref. [[]]. It has been evaluated by varying in the SU(3) chiral expansion the value of the LO low-energy constant 
/ in the range 100-131 MeV (ETMC finds that the effect of a similar change in its chiral fit is smaller than the error 
already assigned to the chiral extrapolation ^5|). Since the updated result (3^) is in very good agreement with the 
lattice average given in eq. (p. a), I find unnecessary to change this average. 
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represent the lattice average (fT8|). The lattice-independent estimate (3.5) of /+(0), obtained using 
the unitarity determination of V us , is also shown for comparison. In the lower side of the plot, we 
collect the quark model result of Leutwyler and Roos (LR 84) [[7|] together with the results of more 



recent analytical model calculations []260-[[29|]. These latter determinations of /+(0) turn out to be 



larger than both the value predicted by lattice QCD and the result implied by CKM unitarity. 

3.2 Lep tonic kaon decays: fk/fn 

A detailed review of lattice results for the pion and kaon decay constants has been given by 



Enno Scholz in his plenary talk at this conference [22]. For this reason, in this talk I will only 
summarize the main features of the new calculations and derive a lattice average for the ratio 
fxlfni which is the relevant hadronic input parameter for the determination of the Cabibbo angle. 

A list of unquenched lattice results for fx/fn is collected in table |] (based on the FLAG 
work), with the colour code rating in the FLAG style assigned for each calculation. The relevant 
simulation parameters for these calculations are given in table ||. 

A first look at tables |] and |5] is sufficient to reveal the large number of unquenched predictions 
for fx/ fn obtained by the various lattice collaborations in the last one or two years. For the purpose 
of this conference, the results which are new or have been updated in the last year are framed in 
the tables. Moreover, several lattice predictions for fx/ fx in tableware rated without red squares, 
indicating that all sources of systematic uncertainties in these calculations are kept sufficiently well 
under control. 

At this conference, an update of the MILC analysis has been presented by Urs Heller [[34|]. 
This analysis is based on SU(3) (rooted) staggered ChPT and includes, as a new feature, the N 2 LO 
continuum chiral logs. The chiral fits are based on results obtained with the "fine", "super-fine" 
and "ultra-fine" MILC ensembles and are performed in two stages. The first one only includes the 
lowest quark masses, and it is used to determine the LO and NLO low energy constants. Once these 
constants are fixed, then the higher order contributions, namely the complete N 2 LO contribution 
together with N 3 LO and N 4 LO analytic terms, are determined through a global fit over all quark 
masses. 

The RBC-UKQCD analysis of fx/fx, presented at this conference by Bob Mawhinney [p3|], 
is based on SU(2) chiral fits. While previous results by the collaboration were obtained at a sin- 
gle value of the lattice spacing (a -1 = 1.72 GeV), the new analysis includes data from a second 
ensemble with a finer lattice (a -1 = 2.32 GeV), which also includes lighter dynamical quarks. In 
addition, on the coarse lattice more configurations have been generated, by almost doubling the 
statistics from earlier works. 



An accurate prediction for fx/fn nas been also presented by the BMW collaboration [ pl| , |32| ], 
based on their extensive simulation performed at three values of the lattice spacing, large volumes 
and simulated pion masses reaching down about 190 MeV. The chiral extrapolation is performed 
testing three different assumptions for the quark mass dependence: ChPT, either SU(3) or SU(2), 
or a simple polynomial extrapolation. With respect to the preliminary determination presented in a 
poster by Alberto Ramos at this conference, the final result for fx/fit given in [HI 32] and quoted 
in table |] has a reduced systematic uncertainty. 

The JLQCD/TWQCD collaboration has presented a preliminary determination of the pion and 



kaon decay constants [36] based on a simulation with Nf ■ = 2 + 1 dynamical overlap fermions, at 
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Collaboration 



Ref. 



N 



f 



27 



.i 



i 



ALVdW 09 




BMW 09 




RBC/UKQCD 09 


MILC 09b 




MILC 09a 





JLQCD/TWQCD 09 



PACS-CS 08 
HPQCD/UKQCD 07 
RBCAJKQCD 08 
NPLQCD 06 
MILC 04 



2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 



ETMC 09 



ETMC 07 

QCDSF/UKQCD 07 



fx I fx 



C 


• 


• « 


» 1.192(12)(16) 


p 


★ 


★ i 


r 1.192(7)(6) 


C 


• 


★ « 


» 1.225(12)(14) 


A 


★ 


★ i 


f 1.198(2)(+ 8 6 ) 


A 


* 


★ i 


r 1.197(3)(+ 6 3 ) 


C 


• 


■ 1 


1 1.210(12W 


A 


* 


■ 1 


1 1.189(20) 


A 


* 


• i 


r 1.189(2)(7) 


A 


• 


* ■ 


1 1.205(18)(62) 


A 


• 


■ 1 


1 1.218(2)0 


A 


★ 


• « 


» 1.210(4)(13) 


A 


• 




r 1.210(6)(15)(9) 


A 


• 


• ■ 


1 1.227(9)(24) 


C 


• 




» 1.21(3) 



Table 4: Colour code rating in the FLAG style and summary of unquenched lattice results for fa/ fn- The 
new results, obtained in the last year, are framed. 



Collaboration 



Ref. 



yv 



/ 



action 



a/fm (Lm n ) min m^/MeV 



ALVdW 09 




BMW 09 




RBC/UKQCD 09 


MILC 09b 




MILC 09a 





JLQCD/TWQCD 09 



PACS-CS 08 
HPQCD/UKQCD 07 
RBC/UKQCD 08 
NPLQCD 06 
MILC 04 



[0 

H 

n 

[0 



IP 



2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 
2+1 



KSmilc/DWF 
impr. Wilson 
DWF 



MILC 
MILC 



MILC 



KS 

MILC 

Overlap 

clover (NP) 

^oHISQ 
^MILC 

DWF 

KSmilc/DWF 



t^qMILC 
^^MILC 



>0.09 
>0.07 
>0.08 

> 0.045 

> 0.045 
0.10 
0.09 

>0.09 
0.11 
0.13 

> 0.09 



3.8 > 300 

4.0 > 190 

4.0 > 290 

4.0 > 230 

3.8 >230 

2.8 > 340 

2.3 > 160 

3.8 > 250 

4.6 > 330 

3.7 > 290 

3.8 > 250 



ETMC 09 



ETMC 07 

QCDSF/UKQCD 07 



[0 
10 

[0 



max. tmQCD 
max. tmQCD 
clover (NP) 



>0.07 
0.09 
> 0.06 



3.2 > 260 
3.2 > 300 
4.2 > 300 



Table 5: Parameters of the simulations listed in table |j. 
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five different up and down quark masses. While the chosen approach is definitively a benchmark 
for lattice QCD, since the overlap formulation preserves an exact chiral symmetry at finite lattice 
spacing, the simulation performed by JLQCD/TWQCD relies yet on a single value of the lattice 
spacing (a ~ 0.10 fm) and a rather small lattice size ((Lm K ) mm < 3), two features which imply the 
red squares in the FLAG stye rating of table |j. 

Among the other Nf = 2 + 1 determinations of fx/ fn collected in table ^ and not presented at 
this conference, the ALVdW 09 and HPQCD/UKQCD 08 results are also free of red tags and will 
be thus included in the determination of the final lattice average. 

A determination of fx/ fn with Nf = 2 dynamical quarks, in which all sources of systematic 
uncertainties are kept under control (i.e. no red squares in the colour code rating) has been recently 
obtained by the ETM collaboration [|T]], using maximally twisted mass fermions. The comparison 
with the most precise Nf ■ = 2 + 1 calculations suggests that the quenching effect of the strange 
quark in the determination of fx/ fx, is currently negligible, within the reached lattice accuracy. 
Therefore, I will include this result in the determination of the final lattice average (in this case, 
given the large number of independent determinations of fx/ fn which are free of red squares in 
table the inclusion of the ETM 09 result has a marginal impact). I then obtain 



/*//* = 1.196(1)(10) 



(3.11) 



where the first error is statistical, evaluated assuming gaussian statistical uncertainties, and the 
second error is systematic. For the latter, I'm quoting an uncertainty which is of the same size 
of the typical systematic uncertainty estimated for the most precise lattice determinations listed in 
table § 

The average ( |3.1l| ) can be translated into a determination of the kaon decay constant, fx, by 
using f % = 130.41 (3) (20) MeV [14], as determined from the measurement of the leptonic pion 
decay rate and V ll( i from eq. (O). In this way one gets 



fx = 156.0 ±1.3 MeV . 



(3.12) 



The lattice average (3. 1 1 ) is also in very good agreement with the lattice-independent estimate 



of fx/ fn quoted in eq. (3.5), based on CKM unitarity and the determination of V u d from nuclear 



j3-decays. By combining the lattice prediction (p.l 1[) with the experimental result of eq. (EO) and 



with Vud from nuclear /3-decays, one obtains the estimate 



1^1^=0.2249(19) . 



(3.13) 



This in good agreement with both the determination (3.10) from Kg} decays and with the unitarity 



determination in eq. ( |3.4| ). 

A summary of lattice results for fx/ fn is shown in fig. [I], which is an updated version for 
the lattice conference of a plot produced by FLAG [[]]]. The blue square and grey band in the plot 



represents the lattice average (3.11). The lattice-independent estimate (3.5) of fx/ fn is also shown 
for comparison. 



13 



Kaon physics from lattice QCD 



Vittorio Lubicz 



L/f 

K 71 

1.14 1.16 1.18 1.20 


1.22 1.24 1.26 


I 

I — I — 
H 

h 

I 

h 


■ h 

-»-H 

I — HH 
— IBH 
I — 

■ 

■ 


I I 

— I ALVdW 09 
BMW 09 

H ■ H RBC/UKQCD 09 

MILC 09b 
MILC 09a 

* 1 JLQCD/TWQCD 09 

H PACS-CS 08 

HPQCD/UKQCD 08 
1 1 1 RBC/UKQCD 07 


I 


1 

1 


■ — 1 NPLQCD 06 

W 1 MILC 04 


h- 


H 


CH 1 ETMC 09 

1 — □ — 1 1 ETMC 07 

O 1 QCDSF/UKQCD 07 




■ 


lattice average 


I 


HH 


nuclear p decay 



1.14 1.16 1.18 1.20 1.22 1.24 1.26 



Figure 2: Summary of unquenched lattice results for the ratio of decay constants fx/ fx as obtained from 
simulations with Nf = 2 and Nf = 3 dynamical quarks (red squares). The black circle represents the lattice- 



independent estimate of fx / f% given in eq. (3.5), which is valid only in the Standard Model. The blue square 



and the grey band show the average of the lattice results for / f„ derived in eq. ( 3.11 ). This plot has been 
produced by FLAG and updated for the lattice conference. 

4. K° - K° mixing: B K 

The accuracy in the lattice determination of the kaon bag parameter, Bk , has remarkably im- 
proved over the last few years. This progress is well illustrated by the following sample of lattice 
averages presented at the lattice conferences: 

Lattice '96 B K = 0.90 ±0.03 ±0.15 S. Sharpe, @ (4.1) 



Lattice '00 B K = 0.86 ±0.06 ±0.14 L. Lellouch, @ (4.2) 



Lattice '05 B K = 0.79 ±0.04 ±0.08 C. Dawson, @7| (4.3) 



Lattice '08 B K = 0.723 ± 0.037 L. Lellouch, Q , (4.4) 

where B K is the renormalization group invariant definition (for Nf = 3) of the bag parameter. The 
second error quoted in eqs. (p~T|)-(p3|) was an estimate of the quenched error. This uncertainty has 
limited the accuracy of the lattice calculations of Bk for a long time. It has started to decrease when 
the first unquenched estimates have been performed in the last few years, and in the average of Bk 
quoted at the last two lattice conferences it was definitively removed. Nevertheless, all unquenched 
determinations of Bk available until last year were obtained at a fixed (and rather large) lattice spac- 
ing. Thus, a quantitative estimate of discretization effects affecting these calculations, which could 
have been not negligible according to the experience gathered in the quenched approximation, was 
not available yet. 

The list of unquenched lattice results for Bk available today is collected in table ^, with the 
colour code rating in the FLAG style assigned for each calculation. The relevant simulation param- 
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Collaboration 



Ref. 



9 / 



/ 



9 



J 



4 



/ 



< s (2GeV) ^ 



ALVdW 09 




RBC/UKQCD 09 


SBW09 





RBC/UKQCD 07 g ||] 
HPQCD/UKQCD 06 pi 



2+1 


A 


• 


★ 


• 


★ 


• 


0.527(6)(20) 


0.724(8)(28) 


2+1 


C 


• 


• 


* 


★ 


• 


0.537(6)(18) 


0.738(8)(25) 


2+1 


C 


★ 


★ 


■ 


■ 


• 


0.512(14)(34) 


0.701(19)(47) 


2+1 


A 


■ 


• 


* 


★ 


• 


0.524(10)(28) 


0.720(13)(37) 


2+1 


A 


■ 


• 


* 


■ 


• 


0.618(18X135) 


0.83(18) 



ETMC 09 1 [j5Vp 2 C ★ • • * • 0.518(21)(21) 0.730(30)(30) 

JLQCD08 HI 2 A ■ • ■ ★ • 0.537(4)(40) 0.758(6)(71) 

RBC04 Q 2 A ■ ■ * • 0.495(18) 0.699(25) 

UKQCD04 Q 2 A • 0.49(13) 0.69(18) 



Table 6: Colour code rating in the FLAG style and summary of unquenched lattice results for Bk- The 
symbol ■' means that these results have been obtained at [M K L) m i n > 4 in a lattice box with a spatial 
extension L < 2 fm. The new results obtained in the last year are framed. 



Collaboration 



Ref. 



Nt 



action 



a/fm {Lm 7! ) min m^/MeV Ren. 



ALVdW 09 



RBC/UKQCD 09 


SBW09 





RBC/UKQCD 07 
HPQCD/UKQCD 06 



2+1 


KSmilc/DWF 


>0.09 


3.5 


>230 


Rl 


2+1 


DWF 


>0.08 


4.0 


>290 


Rl 


2+1 


Ty-cHYP 
^MILC 


>0.06 


2.5 


>200 


VTU 


2+1 


DWF 


0.11 


4.6 


>330 


Rl 


2+1 


K-oMTLC 
a MTLC 


0.12 


4.6 


^360 


PT1£ 



ETMC 09 



JLQCD 08 
RBC04 
UKQCD 04 



2 


tmQCD/OS 


>0.07 


3.2 


^260 


Rl 


2 


overlap 


0.12 


2.8 


>290 


Rl 


2 


DWF 


0.12 


4.6 


>490 


Rl 


2 


clover(NP) 


0.10 


6.2 


>740 


VTU 



Table 7: Parameters of the simulations listed in table I 



eters for these calculations are given in table ^. 

The novelty this year are a number of new calculations, namely ALVdW 09 [49], RBC/ 
UKQCD 09 g, SBW 09 and ETMC 09 [§7|], which are framed in table | for better 

illustration. At variance with the previous unquenched calculations, they all involve an extrapola- 
tion to the continuum limit for B K , based on two (ALVdW 09, RBC/UKQCD 09) or three (SBW 09, 
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ETMC 09) values of the lattice spacing. 

While the RBC/UKQCD 09, SBW 09 and ETMC 09 results are still preliminary, having been 
only presented in the proceedings of this conference, the ALVdW 09 calculation is already pub- 
lished @. It uses a mixed action setup, with domain wall fermions valence quarks over the 
Nf = 2+1 staggered gauge field configurations produced by MILC ("fine" and "coarse" lattices). 
The choice of domain wall fermions for the valence allows a straightforward implementation of 
the RI-MOM method in order to non-perturbatively renormalize the four-fermion operator relevant 
for Bk (in the previous determination of Bk with staggered fermions by HPQCD/UKQCD [p6|], 
the one-loop perturbative determination of the renormalization constant turned out to be affected 
by an uncomfortably large systematic uncertainty). With domain wall fermions, due to the residual 
breaking of chiral symmetry which is allowed by the finite extension of the lattice in the fifth di- 
mension, the Bk operator has a small mixing with operators of wrong chirality (but not with those 
of incorrect taste). The final result for the bag parameter is obtained after performing a combined 
chiral and continuum extrapolation based on NLO SU(3) "mixed action" ChPT t plj |62| ], with the 
inclusion of some analytic N 2 LO contribution. The final accuracy quoted for Bk is about 4%, 
where the dominant source of uncertainty is represented by the determination of the renormaliza- 
tion constant. 

The new result for Bk obtained by RBC/UKQCD, which updates the first precise unquenched 
calculation of ref. [55, 20], has been presented by Chris Kelly at this conference [50]. The main 
improvement with respect to the previous calculation is the use of the second finer lattice ensem- 
ble simulated by RBC/UKQCD, which allows to perform a combined SU(2)-chiral and continuum 
extrapolation of Bk- This, in turn, has permitted a substantial reduction of the systematic error 
( 3.5%), since in the previous calculation discretization effects represented the main source of un- 
certainty. The new calculation also makes use of reweighting in the strange sea sector, with a cor- 
responding interpolation in the valence sector, to reach the physical strange quark mass. Renormal- 
ization of the four-fermion operator is performed non-perturbatively using the RI-MOM approach 
generalized to various non-exceptional momentum renormalization conditions (for a detailed dis- 
cussion of this approach see the plenary talk by Yasumichi Aoki at this conference [63]). 

The new SBW 09 calculation [fl]H|4|] (where the acronym indicates the Seoul, Brookhaven 
and Washington institutions) uses a mixed action setup, with HYP-smeared staggered valence 
fermions, which are effective at reducing taste-breaking effects, and the asqtad staggered sea 
quarks, i.e. the MILC ensembles. Results are obtained at three values of the lattice spacing, 
namely the "coarse", "fine" and "super-fine" MILC ensembles. The chiral extrapolation is per- 
formed by using the proper either SU(3) or SU(2) mixed action staggered ChPT at NLO, with or 
without adding an analytical N 2 LO term. The SU(2) result, which is the one quoted in table ^, leads 
eventually to a smaller error. Since the minimum value of m K L in the simulation is about 2.5, this 
calculation gets a red square in table |6| for finite volume errors. These errors, however, have been 
explicitly investigated by the collaboration and they are found to be subdominant for Bk, smaller 
than 1% for the SU(2) analysis. The main uncertainty in the calculation comes from the evaluation 
of the renormalization constant, which is performed in one-loop perturbation theory. With respect 
to the older staggered calculation of ref. [|56|], this uncertainty is now significantly reduced, mainly 
because of the use of the "fine" and "super-fine" MILC ensembles for which the error coming from 
the truncation of the perturbative series, of order a s (\/a) , is expected to be smaller. Nevertheless 
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it still dominates the systematic uncertainty. The collaboration plans to reduce this uncertainty in 
various way. One is to work on a yet finer lattice, another is to use two-loop matching and, fi- 
nally, by using the RI-MOM non-perturbative method. This latter approach, which looks the most 
promising, has been already successfully applied with staggered fermions to the renormalization 
of bilinear quark operators, as discussed by Andrew Lytle at this conference [|64|]. 

The fourth new result for Bk has been presented at this conference by the ETM collabora- 



tion [p7p. This calculation also uses a mixed action setup, with Nf = 2 maximally twisted sea 
quarks and a suitable Osterwalder-Seiler variant of the twisted mass action for the valence quarks. 
This setup simultaneously guarantees the absence of mixing with operators of wrong chirality, 
which is usually present with Wilson-like fermions, and automatic ^(a)-improvement of Bk [p^]. 
The calculation is performed at three values of the lattice spacing (a ~ 0.07, 0.09, 0.10 fm) and 
the physical value of Bk is eventually reached through a combined SU(2)-chiral and continuum ex- 
trapolation. Renormalization of the four fermion operator is performed non-perturbatively with the 
RI-MOM method. Discretization effects in the evaluation of the renormalization constant, which 
starts at ff(g 2 a 2 ), are further reduced by subtracting from the RI-MOM four-fermion Green func- 
tion the leading contribution, explicitly evaluated in &(a 2 ) one-loop lattice perturbation theory, as 
illustrated by Fotos Stylianou in a poster at this conference [ |66|] . 

In order to derive the final lattice average for Bk, I consider the three results which are free of 
red squares in table ||, namely: 

B K = 0.724(8) (28) [ N f = 2 + 1 , ALVdW 09 ] 

B K = 0.738(8) (25) [Nf = 2+1 , RBC/UKQCD 09 ] (4.5) 
B K = 0.730(30) (30) [N f = 2, ETMC 09 ] . 

The agreement among the above determinations, as well as with the new SBW 09 determina- 
tion [|5l|]-[|54|], is remarkable. It is also worth noticing that the result obtained with Nf = 2 dy- 
namical quarks lies in between the two Nf = 2 + 1 determinations, showing that also for Bk the 
effect of quenching the strange quark is not visible, within the accuracy currently reached by lattice 



calculations. For this reason, I average together the three results in eq. ( |4.5| ) and quote as the best 
lattice estimate of Bk the value 

(4.6) 



£k = 0.731(7)(35) 



While the statistical error is evaluated assuming gaussian statistical uncertainties, the systematic 



uncertainty quoted in eq. (4.6) is slightly increased with respect to the one quoted by the individual 
calculations, to account for the preliminary status of two out of the three calculations on which 
the final average of Bk is based. A summary of the unquenched results for Bk is also presented 
in fig. |3| which is an updated version for this conference of a plot produced by FLAG [ffj. The 



average value of Bk, given in eq. (4.6), is shown in the plot with a blue square and a grey band. 

The Sk parameter which controls the amount of indirect CP violation in K° — K° mixing, 
and whose theoretical determination relies on Bk, plays a relevant role in the unitarity triangle 
analysis, both within and beyond the Standard Model. Since in the Standard Model the analysis is 
largely overconstrained, it can be also used to extract the values of the relevant hadronic parameters, 



including Bk [67]. The latest determination obtained in this way by the UTfit collaboration, which 
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Figure 3: Summary of unquenched lattice results for the RGI parameter Bk as obtained from simulations 
with Nf = 2 and Nf = 3 dynamical quarks (red squares ). The blue square and the grey band show the average 
of the lattice results derived in eq. (4.6). This plot has been produced by FLAG and updated for the lattice 
conference. 



is therefore only valid in the Standard Model, is 

Mm = 0-87(8). (4.7) 

This estimate also takes into account the long distance contributions to both the absorptive and the 
dispersive part of the AS = 2 effective Hamiltonian evaluated in refs. |6^] and summarized in 
the multiplicative factor k e = 0.94 ±0.02. The Standard Model prediction of Bk from the unitarity 
triangle analysis, given in eq. (|47/j), shows a tension with the direct lattice determination (1.6), at 



the level of 1.6 a. Whether such a tension should persist, with the increasing precision of both 
lattice calculations and of the unitarity triangle analysis, it could become a clear signal of physics 
beyond the Standard Model. 

In new physics models, like for instance the supersymmetric extensions of the Standard Model, 
the theoretical expression of e k depends in general on the complete basis of eight AS = 2 four- 
fermion operators [J70|]. Due to parity conservation in the strong interactions, only five of these 
operators have independent matrix elements. The knowledge of these matrix elements is then cru- 
cial in order to derive reliable predictions for Ek in the context of specific new physics models. 
Lattice calculations for the complete basis of AS = 2 four-fermion operators have been only per- 



formed so far in the quenched approximation [71, 72, 73], and the results turn out to be in poor 



agreement among each other. Very preliminary results for the full basis of AS = 2 matrix elements 



have been presented at this conference by ETMC [57]. Given their high phenomenological interest 



it would be helpful if such a calculation were also address by other collaborations using different 
lattice approaches. 
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5. Non leptonic kaon decays 

While not much time (and space) is left to discuss the lattice studies of non leptonic kaon 
decays, I cannot conclude this talk without at least mentioning the extraordinary effort which is 
being put forward by the RBC/UKQCD collaboration in addressing this issue. The results obtained 
last year by the collaboration have shown that the attractive approach which uses ChPT to relate the 
K — > %% matrix element of interest to the simpler matrix elements of K — > % and K — > transitions 
is affected, in the kaon mass region, by significantly large chiral corrections. This is an intrinsic 
uncertainty, which cannot be avoided. Therefore, the direct calculation of K — > nn matrix elements 
on the lattice must be addressed. 

At this conference, three talks have been dedicated to this topic by RBC/UKQCD 0, 



The main indication is that, while a substantial computational efforts will be required in order to 
obtain a reasonably accurate estimate (i.e. at the level of 10-20%), the direct calculation of both the 
Aq and A 2 complex amplitudes is however feasible. The required theoretical and numerical tools 
include, in particular, the use of chiral fermions, non-perturbative RI-MOM renormalization and 
finite volume methods. 

An exhaustive description of this topic would require a dedicated talk by itself. For that, I 
would like to refer the reader to the excellent review given by Norman Christ at the KAON'09 
conference [|77|]. 



6. Summary and outlook 

The number of large unquenched lattice simulations which are being applied to the study kaon 
physics is rapidly increasing. This effort has allowed to achieve in the lattice determination of some 
(relatively simple) kaon physics observables an unprecedented accuracy. In this talk, I concentrated 
most of the attention in reviewing the lattice results for three of these parameters, namely the vector 
semileptonic form factor /+(0), the ratio of decay constants fy/fn and the kaon bag parameter Bk- 
Lattice averages for these parameters have been given in eqs. ([T8|), ( J3.ll ) and (4.6). Lattice studies 



of non leptonic kaon decays are significantly more challenging. Nevertheless, important progress 
has been achieved also in this field, and first, reliable results for both the AI = 1/2 rule and e'/e 
are expected to produced in two ore three years. 
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